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We present the microscopic theory for the BEC-BCS crossover of an atomic Fermi gas in an optical 
lattice, showing that the Feshbach resonance underlying the crossover in principle induces strong 
multiband effects. Nevertheless, the BEC-BCS crossover itself can be described by a single-band 
model since it occurs at magnetic fields that are relatively far away from the Feshbach resonance. A 
criterion is proposed for the latter, which is obeyed by most known Feshbach resonances in ultracold 
atomic gases. 
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Introduction. — Ever since the compelling work of De- 
Marco and Jin 0], the study of degenerate Fermi gases 
has been at the forefront of ultracold atomic physics. 
Much of the impetus behind current research in this field 
was provided by the successful experimental investiga- 
tion in the last two years of the crossover between the 
Bose-Einstein condensate (BEC) of molecules and the 
Bardeen-Cooper-Schrieffer (^CS) state of Bose-Einstein 
condensed Cooper pairs UHli, SS HH- The pair- 
ing observed in an unequal spin mixture is presently un- 
der intense scrutiny [9|, |lOj . Another exciting direction 
being explored lately is the physics of ultracold Fermi 
gases in an optical lattice [ll|. The latter gives, for in- 
stance, the possibility to experimentally solve the famous 
positive-f/ Hubbard model with repulsive on-site inter- 
actions which is believed to embody the physics of high- 
temperature superconductors [l2, UM ■ The negative-f/ 
Hubbard model with an attractive on-site interaction [lj| 
can also be explored. This model is interesting for exam- 
ple because, at half filling, particle-hole symmetry leads 
to an exact SO (A) symmetry. Therefore, the physics at 
half filling is analogous with that of the SO{5) theory pro- 
posed for high-temperature superconductors [isl IigL Il7l| . 
Moreover, at low filling fractions, a BEC-BCS crossover 
takes place. This is the subject of the present Letter. 

In the atomic gases of interest, studying the BEC- 
BCS crossover demands the exploitation of a Feshbach 
resonance [l^ |l9j- As a result, the physics is much 
richer than in the negative-t/ Hubbard model and re- 
quires a multiband description. This is elucidated in 
Fig. n In the absence of an optical lattice, and includ- 
ing solely two-body effects, a dressed molecule formally 
exists only below the Feshbach resonance with a binding 
energy which, near resonance, depends quadratically on 
the magnetic field. Its binding energy is shown as the 
dash-dotted line. Above the Feshbach resonance, the 
dressed molecule is merely a resonance state as it then 
has a finite lifetime due to the possibility of decaying into 
the atomic continuum [20j as represented by the dotted 
line continuation. 

Placing the system in a sufficiently deep optical lattice. 
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FIG. 1: Comparison of the dressed molecular binding energy 
at various magnetic fields for the homogenous, the on-site 
two-body, and the many-body cases corresponding to the 
dash-dotted, dashed, and solid lines, respectively. We use 
K atoms near the Feshbach resonance at Bo = 202. IG and 
with a filling fraction of 0.1. The on-site harmonic oscil- 
lator frequency a; is 27r x 58275Hz, which corresponds to a 
lattice with wavelength A = 806nm and a Rabi frequency of 
fin = 2n X 1.43GHz. In the homogeneous case, the molecules 
become short lived resonances above Bo and this is indicated 
by a dotted curve. The inset shows the two-particle con- 
tinuum associated with the lowest Bloch band. The two- 
body binding energy crosses through twice the Fermi energy 
at B ~ SOOG, demonstrating that the BEC-BCS crossover 
occurs far from Bo- The upper curve in the inset corresponds 
to the hole-particle symmetric solution of the BCS gap equa- 
tion. 



the bare (closed-channel) molecular state of the Feshbach 
resonance interacts primarily with a discrete set of atomic 
states instead of an atomic continuum. Thus, a dressed 
molecule now exists for all magnetic fields and its en- 
ergy, plotted as dashed lines, shows a number of avoided 
crossings before it ultimately becomes a resonance in the 
continuum above the optical lattice potential. 

Using this result, we can restate the interacting prob- 



lem for deep optical lattices in terms of atoms and dressed 
molecules emerging from the on-site two-body solution. 
The binding energy of these dressed molecules obtained 
from our exact solution of the on-site Feshbach prob- 
lem is in good agreement with recent experimental data 
|2lL 1221 , showing the validity of this approach. Many- 
body physics enters if we allow both the atoms and the 
dressed molecules to hop to adjacent lattice sites |2,Sj . 
In particular, for the low filling fractions considered, the 
two-body binding energy of the dressed molecules passes 
through the Fermi sea in the lowest Bloch band, as shown 
in the inset of Fig. ^ It is due to this phenomenon that 
the BEC-BCS crossover occurs in this resonantly inter- 
acting case [2^. The interaction of the molecules with 
the Fermi sea further clads the dressed molecules with 
many-body effects resulting in a change of the binding 
energy such that it only asymptotically reaches twice the 
Fermi energy at high magnetic fields, as shown by the 
solid line. Mathematically, this arises because the two- 
body pole at 3fku in the self-energy of the molecules is 
replaced by the famous logarithmic BCS singularity at 
twice the Fermi level. By placing the system in an op- 
tical lattice, we thus find that the BEC-BCS crossover 
takes place rather far from the resonance. Consequently, 
as we shall show below, we may always use the single- 
band approximation to describe the BEC-BCS crossover 
for all experimentally relevant Feshbach resonances. 

BEC-BCS crossover theory in the lattice. — We now 
begin with an outlay of our general BEC-BCS crossover 
theory in the lattice before discussing the specific case 
of ^"K. The second-quantized grand-canonical Hamilto- 
nian in the presence of an optical lattice potential V^(x) 
is 
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where V'o-(x) are the atomic and V'm(x) the molecu- 
lar annihilation operators, and fj, is the chemical po- 
tential. The atom-molecule coupling is given hy g — 
?i-\/47rabgAi?A/iniag/"T^, where abg is the background 
scattering length, AB is the width of the Feshbach res- 
onance, A/imag is the difference in magnetic moments 
between the closed and open channels of the Feshbach 
resonance, and m is the atomic mass [20|, |2J|. The de- 
tuning from resonance is (5 = Afiamg^B — Bq), where Bq 
is the value of the magnetic field B at the resonance. It 



gives the location of the bare molecular level with respect 
to the atomic continuum in the absence of the lattice. 

With a lattice potential present, we turn to a descrip- 
tion in terms of Bloch bands since plane wave states no 
longer diagonalize the Hamiltonian. Summing the effects 
of nearest-neighbour tunneling leads to a dispersion in 
the lowest band of the optical lattice equal to £(k) = 
~2t^[cos{k^X/2)+cos{kyX/2)+cos{k,X/2)]+3fiLj/2,with 
A the wavelength of the lattice laser and 3fiuj/2 the 
on-site, ground-state energy of a single atom. The 
atomic and molecular tunneling amplitudes are given by 

ia,m = ^Vo^Ef^Jir^/'eM-^^n/E^), where ^0 is 
the peak-to-trough depth of the optical lattice potential 
and the atomic and molecular recoil energies are E^ — 
2{-Knf/m\^ and E^ = Ef/2. Moreover, the effective 
on-site atom-molecule coupling g^ — g^ J(ix|xfc(x)|^ 
becomes band-dependent, where Xfc(x) is the Wannier 
function in band b = 0,1,2,... of the optical lattice. 
In the lowest band, we write this effective coupling as 
g' = go — g/{2'KpY^'^ where I = yjhjmio is the on-site 
harmonic oscillator length [23|. 

Obtaining pertinent thermodynamic quantities for this 
system involves the calculation of the molecular self- 
energy. We obtain this self-energy by first integrat- 
ing the fermions out of the partition function exactly. 
We then perform an RPA or Bogoliubov approxima- 
tion around the mean field A = g'iipm) = ff'v^Sc 
of the bare molecular condensate. In the on-site two- 
body limit, we have ti^ — and A = 0, leaving dis- 
persionless harmonic oscillator bands with energies Sb = 
(6 + 3/2)fnjj. The two-body molecular self-energy then 
becomes fiS{E) = g'^y/TrG{E ~ 3fiu!)/fiuj where the an- 
alytic sum over dispcrsionless bands is encapsulated in 
the function ^] G{E) = T{-E/2huj)/T{-E/2liLu-l/2). 
The binding energy of the dressed molecules is then found 
from E = S + 3huj/2 + nE{E). 

In the many-body case, where ia, ^m, and A are 
all nonzero in general, we work in a regime where the 
temperature is sufficiently below the Fermi tempera- 
ture that we may use a zero-temperature approxima- 
tion. Due to the large band gap in a deep optical 
lattice, the structure of all bands but the lowest has 
a negligible effect and they can be taken to be dis- 
pcrsionless, as in the two-body case. However, they 
can not be excluded altogether from calculations since, 
as we have just witnessed, they renormalize the two- 
body energy levels. The BCS dispersion relation of 
the fermions is now fiw(k) — ^(£(k) — fi)^ + |Ap and 
we can write the normal and anomalous molecular self- 
energies in terms of the usual coherence factors of BCS 
theory |u(k)|2 = [fitj(k)+e(k)-^]/2fiw(k) and |i;(k)|2 = 
[fitj(k) - (e(k) - ^)]/2fcj(k) to obtain 
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where k± = k' ± k/2, i? = ifcj„ + 2/i — Sfej, A^g is the 
total number of lattice sites and the sums are over the 
first Brillouin zone only: k'^,k'y,k'^ £ (— 27r/A, 27r/A). An 
equation for the BCS gap A is then found in the form of 
a Hugenholtz-Pines relation 
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fiSn(O,0)-fiSi2(O,0). (4) 



Having these quantities in hand enables us to de- 
termine the partition function of the system from the 
Gaussian effective action of our RPA theory and, via 
the resulting thermodynamic potential 51, we can cal- 
culate the total number of particles using the identity 
N = —dfi/dfi. This leads to the equation of state for 
the total atomic filling fraction n = N/Ng given by 



Tr [rgG^ 



g'- 



Tr [G^] + i Tr 



Q„ 



9/1 



(5) 

where T3 is the third Pauli matrix and the 2x2 Nambu 
space matrices G^^g , C^ and S are the BCS atomic 
Greens function, the molecular Greens function, and self- 
energy matrix composed of Eqs. Q and © respectively. 

Results and discussion. — Although our approach is 
quite general, the numerical calculations in this Letter 
are for the Bq = 202. IG Feshbach resonance of ""^K 
atoms, which has a width of AB — 7.8G [3|, and for 
a total filling fraction of 0.1 as an illustrative and ex- 
perimentally relevant case. For every S, solving the BCS 
gap equation, Eq. Q , and the equation of state simulta- 
neously yields self-consistent values for fi and the filling 
fraction of the bare molecular condensate nj^^ = A^/g'^. 
Since fiuctuation effects are expected to be less important 
away from resonance, in the first instance we neglect the 
fluctuation contribution reflected by the last two terms 
on the right-hand side of Eq. (O. In order to examine 
nj^j,, it proved requisite to include the structure of the 
higher bands to which intermediate states strongly couple 
in the neighbourhood of the Feshbach resonance. Having 
done so, the many-body result monotonically approaches 
the 2-body limit in the resonant region as more bands 
are included in the self-energy calculation. This is clear 
from Fig. 13 When including higher bands, we used a 
multiband generalization of Eqs. ^ and © that did not 
factor in any interband coupling other than the already 
present interaction of the mean field Ab = gb{'>Pni) with 



the molecular condensate. Although this leads to only 
qualitative results in the higher bands, the correct two- 
body physics is obtained and the importance of multi- 
band effects near resonance is well demonstrated. 

We also determine the filling fraction of Bose-Einstein 
condensed dressed molecules as this will provide us with 
another good indication of where the crossover from a 
Bose-Einstein condensate of dressed molecules to the 
BCS state composed of Bose-Einstein condensed Cooper 



pairs takes place. This fraction is given by rii, 



JZ, 



where Z is the so-called wavefunction renormalisation 
factor which determines the probability amplitude for the 
dressed molecular wavefunction to be in the bare molec- 
ular state |2i|. This is plotted in the inset of Fig.Ol For 
this purpose, we did not calculate the true many-body 
Z, which would involve the full fluctuation calculation 
|26j, but used instead the 2-body wavefunction renor- 
malisation factor given by Z^^ = 1/(1 - dHY,{E)/dE), 
where hT,{E) is the two-body self-energy discussed ear- 
lier. This results in an upper bound for n^xc, owing to the 
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FIG. 2: The filling fraction of the bare molecular conden- 
sate, 2n||(./n, drawn as solid lines, decreases monotonically to 
the two-body limit as progressively more bands are included, 
from 1 to 300 in this case. The dashed line is Z^^ , which is 
equivalent to the filling fraction of the molecular condensate 
in the BEG region shown here. 

fact that Z^^ is always smaller than its many-body coun- 
terpart Z . The magnetic field region where this dressed 
filling fraction decreases from one to zero indicates the 



position of the BEC-BCS crossover, which is clearly in 
the same neighbourhood as where the two-body molecu- 
lar energy enters the Fermi sea. From the inset in Fig.O 
we see that the BEC-BCS crossover is dominated by con- 
tributions from the lowest band only; including higher 
bands in the self-energy calculation does not bring about 
any apparent change. 

In Fig. 13 we compare Z^^ with 2n^^/n. Indeed, 
2nJ^j./n and Z^^ are similar in the BEC limit, since the 
gas then consists solely of a Bose-Einstein condensate 
of dressed molecules. This is a result of the fact that 
in this regime the energies of the atomic states lie far 
above the dressed molecular state, as seen in Fig.^ which 
suppresses many-body effects. The BCS region, how- 
ever, sees a pronounced difference between the two-body 
renormalisation factors and the bare molecular conden- 
sate fraction. 
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FIG. 3: Two-body wavefunction renormalisation factor Z^^ 
(dashed line) plotted together with the filling fraction of the 
bare molecular condensate 2n^^/n including a varying num- 
ber of bands (solid lines), which shows that a single-band 
calculation is adequate in the BEC-BCS crossover region of 
interest. The inset shows the approximate filling fraction of 
dressed molecules. Time, in the crossover region. The multiple 
solid lines of the main graph all lie on top of each other in 
this high magnetic field region. 

Our solutions of the BCS gap equation confirm the 
presence of the logarithmic BCS singularity arising at 
high detuning where the gap then scales as A ~ 
2t'"exp(-29t°-/|t/cff|) Hi. Here, U^s ^ - g'^ / (S + 3huj - 
2/i), is the on-site Hubbard attraction of the atoms. In 
the crossover region, the molecular energy crosses the 
lowest-band boundary when S ~ g'^ /\2ts,. Furthermore, 
C/cff indicates this crossing occurs close to a Feshbach 
resonance when 5 ~ ifiuj/2. This leads to the criterion 
AS ^ 18ia \/47rfi/77T,u;/abg A//inag which ensures that the 
Feshbach resonance and the BEC-BCS crossover are well 
separated, as required for the applicability of a single- 
band model to the BEC-BCS crossover in an optical lat- 



tice. Most known Feshbach resonances adhere to this 
criterion. 

In conclusion, we have presented the full micro- 
scopic many-body theory that describes the BEC-BCS 
crossover in an optical lattice due to the presence of a 
Feshbach resonance. Both two-channel physics and the 
possible effects of higher bands were taken into account. 
The many-body Z can only calculated from the full fluc- 
tuation theory and this is an important topic of further 
study. 
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